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CALCULATION  OF  THE  POWER  DENSITY  IN  THE 
FREDHOLM  EQUATION  THAT  YIELDS  A 
WEIBULL  pdf 


1.  INTRODUCTION 

Detection  in  non-Gaussian  noise  is  a  very  difficult  subject  to  study  in  radar.  To  maximize  detec¬ 
tion  probability,  the  optimum  filter  is  the  Neyman-Pearson  likelihood  ratio  test.  It  is  a  difficult  prob¬ 
lem  to  obtain  the  multivariate  probability  density  function  required  foi  generating  the  optimum  niter. 
Several  detectors  for  detecting  a  signal  in  non-Gaussian  noise  have  been  found  by  Cantrell  [  1 ,2]  and 
by  Martinez  et  al.  [3].  These  detectors  were  obtained  by  applying  the  Neyman-Pearson  test  after  a 
closed-form  solution  of  the  multivariate  probability  density  function  had  been  found. 

This  report  considers  clutter  that  has  a  non-Gaussian  distribution.  Specifically,  clutter  has  a 
greater  probability  of  having  a  large  value  than  the  Rayleigh  model  obtained  from  a  Gaussian  distribu¬ 
tion,  and  thus  it  requires  a  higher  tailed  distribution  to  model  it  properly.  This  higher  tail  must  be 
properly  accounted  for  in  the  design  of  a  detector  to  avoid  false  alarms.  Generally,  the  Rayleigh  pdf 
underestimates  the  amplitudes  obtained  from  real  clutter.  One  physical  model  that  has  been  proposed 
to  describe  the  non-Rayleigh  nature  of  sea  clutter  amplitude  statistics  is  the  composite  surface  scatter¬ 
ing  model  [4-6]  that  describes  the  fluctuation  of  clutter  amplitudes  by  a  conditional  Rayleigh  pdf,  con¬ 
ditioned  on  a  varying  clutter  power  level.  The  overall  non-Rayleigh  pdf  for  the  clutter  amplitude  is 
given  by 


g(l*  I)  =  j0  exP 
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where  g  ( |  x  | )  is  the  non-Rayleigh  pdf  of  clutter  amplitude  return  |  x  |  and  /  (a2)  is  the  pdf  of  clutter 
power  level  a2.  This  is  a  Fredholm  equation  of  the  first  kind  for  the  unknown  pdf  /  (a2).  Equation 
(1)  will  be  shown  to  be  in  the  form  of  a  Laplace  transform. 


In  principle,  the  averaged  clutter  pdf  /  (a2)  can  be  found  by  inverting  the  transform  if  g(M) 
is  known.  Selecting  g(\x  |)  to  be  a  Weibull  distribution  is  expected  to  better  represent  real  clutter. 
Inverting  the  transform  when  g(  |  jc  | )  is  a  Weibull  density,  is  very  difficult,  so  an  approximate  solu¬ 
tion  was  found.  The  inverse  was  approximated  by  employing  the  method  of  steepest  descent  to  evalu¬ 
ate  the  Laplace  inverse  contour  integral. 


In  Section  1 ,  we  discuss  the  nature  of  a  spherically  invariant  random  process  (SIRP)  and  derive 
Eq.  (1)  by  noting  that  the  clutter  envelope  is  a  SIRP.  An  asymptotic  expansion  technique  is  applied  to 
the  inversion  of  the  Laplace  transform.  In  Section  2,  we  discuss  the  difference  between  our  results, 
which  is  a  mixture  of  Rayleigh  pdfs,  and  the  Weibull  pdf.  Also  we  generate  random  samples  from 
the  mixture  of  Rayleigh  pdfs  and  the  Weibull  pdf,  and  we  compare  the  result  of  these  two  pdfs. 
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2.  NATURE  OF  A  SIRP  AND  DERIVATION  OF  EQ.  (1) 

A  SIRP  can  be  expressed  as  a  random  mixture  of  n  -dimensional  Gaussian-probable  densities. 
Given  a  random  vector  x  =  (x{  x2  xn )  from  the  random  process,  RP  =  [X (r ),  t(T]  with  mean 
value  n  and  covariance  matrix  K,  a  necessary  and  sufficient  condition  for  the  random  vector  to  be 
from  a  SIRP  is  for  the  probability  density  to  be  of  the  form 

_  IL  - 1 

Pn(x )  =  (2 7T)  2  j  A'  I  2  hn[  (x  -  n)T  K~'  (x  -/!>], 

where  /?„(.)  is  a  quadratic  form.  It  can  be  shown  [6]  that  the  pdf  of  the  clutter  envelope  associated 
with  a  SIRP  must  satisfy 
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We  now  assume  a  Weibull  pdf  for  the  clutter  envelope;  i.e., 
g(|*  |)=  a  in  (2) 
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Substituting  Eq.  (3)  into  Eq.  (2)  and  letting  t  =  —  and  s  =  1  1  yields 
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J  e  sl  h(t)dt  =  A  s  2  exp(-fis2) 


(4) 


This  is  in  the  form  of  a  Laplace  transform.  An  inversion  of  the  Laplace  transform  in  Eq.  (4)  yields 
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a  =  2 


We  can  easily  show  that  the  clutter  amplitude  pdf  Eq.  (3)  is  Rayleigh  for  a  =  2.  In  this  report, 
we  arc  interested  in  the  cases  0  <  a  <  2,  since  we  are  interested  in  a  heavier  tailed  distribution  than 
the  Rayleigh  distribution. 

To  solve  the  inversion  of  the  Laplace  transform  (Eq.  (5)),  we  derived  an  asymptotic  expansion 
of  the  integral  to  find  an  approximation  solution.  Some  sufficient  conditions  that  should  be  met  are 
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discussed  in  the  Appendix.  Since  in  our  problem  all  sufficient  conditions  are  met,  we  can  generate  an 
approximation  solution.  In  Eq.  (5)  for  the  case  0  <  a  <  2  ,  let 


a 


(6) 


Then 


h(t)  =  — — q(s)  exp  {-tp(s))ds.  0  <  a  <  2  (7) 

27 n  Jc-iac 

If  the  singularities  of  <^(i)  all  take  the  form  of  isolated  poles  and/or  isolated  branch  points,  then  by  a 
suitable  deformation  of  the  inversion  contour  we  may  reduce  the  integral  to  the  sum  of  residues  at  the 
poles  plus  a  sum  of  loop  integrals  around  the  branch  points.  A  contour  is  shown  in  Fig.  1.  Equation 
(7)  reduces  to  an  integral  on  the  interval  ( -  oc,  +  oo). 

h(t)  =  — — f  q  (5 )  exp  ( -  tp  (5  ))ds ,  0  <  a  <  2. 

2m  ^  » 
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00 
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The  peak  value  of  the  factor  e~,p(s)  is  located  at  ,v0  the  minimum  of  p(s).  When  t  is  large  this  peak 
is  very  sharp,  and  the  overwhelming  contribution  to  the  integral  comes  from  the  neighborhood  of  sQ. 
It  is  therefore  reasonable  to  approximate  p(s)  and  q(s)  by  the  leading  terms  in  their  power  series 
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expansion  at  s0.  I(t)  is  evaluated  by  extending  the  integration  limits  to  (-00,00)  .  The  result  is 
Laplace’s  approximation  to  /(/).  This  solution  will  also  be  used  for  small  t  and  any  discrepancies  in 
the  resulting  pdf,  and  the  assumed  Weibull  one  will  be  noted. 

We  will  expand  p(s )  in  a  Taylor-series.  The  first  two  terms  will  be  sufficient  since  e~,p{s)  has  a 
sharp  peak  value  at  s0.  One  obtains 


We  will  now  evaluate  the  right-hand  integral.  Let  J=  —  p  (s0)  and  x  =  s  —  s0.  Since 


the  integral  in  Eq.  (9)  equals 


V7>  <■*„) 

ds 


Since  p'(s)  =  0  yields  the  minimum,  from  Eq.  (6)  s0  = 


At  s0  we  have  q(s0)  = 
general  solution  becomes 


.  After  substituting  for  the  value  of  s0  and  p  (5)  in  Eq.  (9),  the 


The  probability  density  function  is  given  by 
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Letting  U  =  - ,  the  probability  density  function  of  clutter  power  that  yields  a  Weibull  density 
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Given  (J  =  — ,  one  more  step  to  simplify  the  equation  becomes 
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3.  DISCUSSION 


In  Section  1,  we  obtained  the  probability  density  function  of  the  clutter  power  level  that  allows 
Wcibull  clutter  to  be  modeled  as  a  SIRP.  The  SIRP  may  be  used  to  describe  the  non-Rayleigh  distri¬ 
bution  including  the  correlation  properties. 

The  approximation  solution  from  Section  1  shows  that  the  Weibull  pdf  may  be  modeled  as  an 
infinite  mixture  of  Rayleigh  pdfs.  To  examine  the  accuracy  of  this  result,  we  compare  the  pdf  gen¬ 
erated  by  numerically  solving  Eq.  (1)  to  the  Weibull  pdf.  Specifically,  we  compare  the  probability  of 
exceeding  |jc  | ,  for  the  two  pdfs,  for  various  values  of  parameter  alpha  between  0  and  2,  any  value 
of  median,  and  any  value  of  jjr  |  . 

The  actual  values  of  Weibull  and  the  Rayleigh  mixture  are  observed  to  be  slightly  different  at 
the  beginning  of  the  curves  as  shown  in  Fig.  2.  We  integrated  the  closed-form  solution  (Eq.  (1)) 
from  zero  to  infinity.  This  integration  was  performed  numerically  by  subdividing  the  integral  into 
small  intervals  and  summing  the  appropriately  weighted  values. 

For  alpha  =  0.5,  the  mixture  of  Rayleigh  pdfs  is  lower  than  the  Weibull  pdf,  since  a  large 
enough  upper  limit  was  not  used  for  the  numerical  integration.  However,  the  curve  can  be  made  to 
approach  the  Weibull  pdf  if  we  use  a  good  upper  limit.  Even  though  this  closed-form  solution  is 
based  on  an  approximation,  it  is  seen  to  be  very  close  to  the  Weibull  pdf,  especially  when  a  =  1, 
where  it  matched  almost  perfectly.  The  pdf  is  particularly  good  at  the  tail  of  the  distribution.  For 
values  of  3  or  greater  of  |jr  | ,  the  two  curves  overlapped  each  other  all  the  way  to  a  probability  of 
I(T7. 


We  also  generated  histograms  based  on  the  two  pdfs.  To  generate  a  histogram,  we  used  a 
transformation  method  to  generate  the  desired  distribution  of  a2,  solve  the  differential  equation,  and 
take  the  inverse  of  the  fu  iiCuOii  .  It  is  necessary  first  to  generate  a  sample  from  the  power  level  pdf 
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/  (<r),  then  to  generate  a  sample  from  the  Rayleigh  pdf  by  using  the  same  a2  of  the  power  level  pdf. 
The  histogram  that  results  from  these  samples  is  the  histogram  of  the  mixture  of  Rayleigh  pdfs. 

Histograms  with  10.000  samples,  are  shown  in  Figs  3  and  4  for  cases  of  a  =  0.75  and  a  =  1.0 
respectively,  along  with  the  Weibull  pdf  histogram  for  10,000  samples.  The  mixture  of  Rayleigh  pdfs 
histogram  is  seen  to  be  very  similar  to  the  Weibull  pdf  histogram.  Also,  we  calculated  the  cumulative 
distribution  function  of  power  level  pdf,  which  indicates  that  the  area  under  the  curve  is  0.999302.  It 
is  very  close  to  I .  In  Figs.  5  to  7,  we  show  the  density  function  for  the  closed-form  solution  and  the 
Weibull  pdf. 

These  results  could  be  used  for  the  problem  of  detecting  targets  in  non-Gaussian  noise,  as 
described  by  Cantrell  11,2],  and  by  Martinez  [3|.  To  obtain  similar  results,  a  closed-form  non- 
Gaussian  i  jltivariate  pdf  is  required.  This  requires  a  closed-form  expression  for  the  density  that  is 
obtained  by  averaging  the  Gaussian  multivariate  over  the  power  density  function  found  in  this  report. 
This  problem  is  still  open  for  solution. 
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4.  SUMMARY 

An  analytical  solution  for  a  mixture  of  Rayleigh  pdfs  of  the  averaged  clutter  return  is  found  y 
using  an  asymptotic  expansion  method.  The  mixture  of  Rayleigh  pdfs  Fits  the  amplitude  fluctuatioi 
of  the  'Weibull  density  It  is  a  very  good  fit  at  the  tail  of  the  distribution,  a  fact  that  allows  a  detec 
tion  threshold  to  be  set  for  a  probability  of  false  alarm  equal  to  10~7. 
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APPENDIX 


The  expansion  technique  described  in  Ref.  16  is  an  effective  way  of  deriving  the  asymptotic 
expansion  of  an  integral  containing  a  parameter.  A  general  type  of  integral  amenable  to  the  method 
of  approximation  is  given  by 

I(x)  =  exp  ( —xp (t ))q (t )dt  ,  (A-l) 

where  x  is  a  positive  parameter  and  the  function  q(t)  is  independent  of  the  parameter  x .  The  peak 
value  of  the  factor  e~xpU)  is  located  at  the  minimum  t0,  and  almost  all  the  contribution  to  the  integral 
comes  from  the  neighborhood  of  r0.  Equation  (8)  is  exactly  the  same  as  Eq.  (A-l).  Thus  it  is  reason¬ 
able  to  approximate  p(t)  and  q  (r )  by  leading  terms  in  their  ascending  power  series  expansions  at  f0. 
Also,  we  need  to  show  that  q(t)  is  infinitely  differentiable  by  taking  the  derivative  in  [  0,oo  )  for  the 
sufficient  conditions.  By  taking  the  derivative  infinitely  from  Eq.  (6),  we  still  get  the  function  of  g 
that  satisfies  one  sufficient  condition. 


q(n\t)  =  0(e*),  0  <  t<  oo 


where  0(eu)  is  an  expression  of  infinitely  differentiable  function.  Also  X  is  a  real  constant  that  is 
independent  of  t .  The  integral  converges  when  x  >  X  and  equals 
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x  0 


n2(x  -  X) 


and  n  is  an  arbitrary  nonnegative  integer.  Therefore 
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I(X)  ~  £ 


;  =0  X 
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since  the  maximum  value  of  |  q{n)(t)  \  is  attained  at  t=  t0 


|cnC*)|  =  |<7<n,('o)l*  ” 

Calculating  the  error  terms  by  using  the  above  equation  yields 

€„(*)  -  en  +  ](x)  =  qin) x'n~\ 
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